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Abstract 

Multi-object astronomical adaptive-optics (MOAO) is now a mature wide-field observation mode to enlarge 
the adaptive-optics-corrected field in a few specific locations over tens of arc-minutes. 

The work-scope provided by open-loop tomography and pupil conjugation is amenable to a spatio-angular 
Linear-Quadratic Gaussian (SA-LQG) formulation aiming to provide enhanced correction across the field 
with improved performance over static reconstruction methods and less stringent computational complexity 
scaling laws. 

Starting from our previous work |lj|, we use stochastic time-progression models coupled to approximate 
sparse measurement operators to outline a suitable SA-LQG formulation capable of delivering near optimal 
correction. Under the spatio-angular framework the wave-fronts are never explicitly estimated in the volume, 
providing considerable computational savings on 10 m-class telescopes and beyond. 

We find that for Raven, a lOm-class MOAO system with two science channels, the SA-LQG improves the 
limiting magnitude by two stellar magnitudes when both Strehl-ratio and Ensquared-energy are used as 
figures of merit. The sky-coverage is therefore improved by a factor of '^5. 

OCIS codes: (000.0000) General. 
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1. Intro: Multi-object adaptive optics systems 


Multi-object Adaptive Optics (MOAO) is now a well-established AO concept to enlarge the 
corrected held to several arcminutes when only a few directions of interest exist in the held 


and need be corrected for 




Similarly to multi-conjugate AO (MCAO) it relies on tomography to provide a 
tridimensional wave-front estimation above the telescope and thereupon provide correction 
on several scientihc relevant directions, reaching out for a multiplexing factor of up to 20 on 
ELT-sized systems 


On a previous paper [l| we have outlined the static minimum-mean square error 
(MMSE) spatio-angular (SA) wave-front estimation (time-dependent wave-front evolution 
not considered) when a Zernike polynomial basis is used and suggested predictive models 
to overcome temporal lag errors when aiming for increased sky-coverage [?]. In the SA 
case estimating the tomographic phase explicitly is circumvented, only the aperture-plane 
is estimated to be later least-squares htted onto the DM - the Learn&Apply approach in 
is similar albeit in a diherent vector space. 

We now move one step forward: we provide extensions to use a zonal basis set coupled 
to the optimal minimum residual phase error variance Linear-Quadratic-Gaussian (LQG) 
controller tailored to MOAO systems.The LQG controller has been presented in the AO 


context in several works 


flls 


lol l and here it is designed and optimized for a MOAO system 


performing open-loop tomography and single-layer optical conjugation. We further propose 
a minimal state representation together with hrst-order near-Markovian predictive models 
which lead to a computational complexity of the order of the static MMSE reconstructors 
for a multiplexing capability of 20 science targets. 

The use of a spatio-angular framework has several advantages over the formulation that 
estimates the layered phase explicitly: 


1. the real-time computation is independent from the number of atmospheric layers; 

2. the analytical spatio-angular covariance matrices can be computed explicitly thus 
avoiding interpolation through the ray-tracing operation to obtain the integrated pupil- 
plane wave-front prohle; 

3. in principle a hybrid data-driven model-based controller can be designed since a subset 
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of the covariance matrices can be acqnired directly from the system (i.e. empirically 
acqnired) instead of relying fnlly on models; alternatively meta-parameters can be 
identihed from which enhanced models are obtained. 


4. the sparseness of the LQG compounding matrices can still, albeit to a lesser degree, 
be exploited for real-time savings as noted in B Q- 


The goal of this paper is therefore two-fold: 

i) fully outline the zonal LQG minium-variance tomographic controller design to optimize 
optical performance, particularly in poor signal-to-noise ratio regimes such that sky-coverage 
can be improved upon by guiding on fainter sources and 

a) provide a suitable. Extremely Large Telescope (ELT)-comphant spatio-angular 
formulation for the Kalman estimator (and static MMSE) that decouples the wave-front 
reconstruction from the number of atmospheric layers used to model it. 

We present results from both simulation and the optical bench of the Raven system, a 
MOAO science and technology demonstrator installed on the 8 m Subaru telescope j^. 

This paper is organized as follows. Section 2 outlines the LQG formulation, section 3 
discusses other SA implementations, namely the static SA reconstructors obtained from the 
LQG, section 4 provides sample numerical results whereas section 5 overviews the Raven 
optical bench and the results obtained in the lab. Section 6 provides the summary and 
conclusions. 


2. Spatio-angular LQG formulation 
2.A. Definitions and assumptions 

We dehne the pupil-plane wave-front indexed by the bi-dimensional spatial 

coordinate vector p = [px, Py] in direction 6 = [6x,0y\ at time t under the hypothesis that 
the turbulent atmosphere is a sum of L thin layers located in a discrete number of different 
altitudes hi as 

L 

'il^{p,e,t) = [li(p]{p,t) = '^uji(pi{p + hie,t) (1) 

1=1 

where (pi{p,t) is the /*^-layer wave-front, uji is the layer relative strength and H is a 
geometric optics propagation operator in the near-held approximation (commonly called 
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ray-tracing 12|]) that relates the aperture-plane wave-front (WF) to the layered wave-fronts 
by adding and interpolating + hiO,t) on the aperture-plane (ho = 0) computational 
grid 131. 

In the following we assume the frozen-flow hypothesis whereby t + r) = <p{p — vr, t) 
with V = [vx'i Vy\ the wind velocity vector. A discrete-time, point-wise representation of the 
WF is also assumed throughout this paper. 

We further assume that averaged phase values over the wave-front sensor (WFS) 
integration period T* are taken at face value, i.e 

1 rkTs 


¥^fc(p) = 


(p{p,t) dt. 


( 2 ) 


J(k-l)Ts 

and that the commands are piece-wise constant over an entire frame (either synchronous or 
not with the WFS sampling frames) 


Ufc = u(t), kTs <t < {k + l)Ts 


( 3 ) 


Ensemble averaging is represented by (•) with shorthand covariance matrix notation 
'^u,v = {uv'^'). Angles are notated in general 0 and guide-star and science directions use ck 
and j3 respectively. 

2.B. LQG regulator synthesis 


The derivation of the LQG for AO applications has been presented in detail elsewhere for 
either the inhnitely fast DM response or otherwise jsl. In this paper we therefore take a 
shorter path for the sake of conciseness and refer the reader to the bibliographic references 
on this particular subject for a more in-depth derivation. 

The discrete-time LQG regulator minimizes the cost function 

M-l 


J(u) = lim — (x^Qx-|-u^Ru-|-2x''^Su) 

^ ^ M^oo 


( 4 ) 


fc =0 


where the weighting triplet {Q,R, S} is made apparent and will be specihed next by 
developing a MOAO-specihc quadratic criterion subject to the state-space model 


Xfc+i = Mxfc -h Buk + Vvk 
Sc,k = Cxk + + Pk 


( 5 ) 
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where is the state vector that contains if not the phase directly a linear combination 
thereof, Sct,k are the noisy measurements provided by the WFS in the GS directions, Uk 
and are spectrally white, Gaussian-distributed state excitation and measurement noise 
respectively. In the following we assume that are known and 

that = 0. 

2.C. Minimum residual phase variance after correction 

The objective cost function considered for MOAO is the minimisation of the aperture- 
plane residual phase variance for individual science directions (3^ G leading to the 

maximisation of the Strehl-ratio. 

If we denote the term l3i,kTs) — l3i^kTs) the residual phase prohle 

after DM correction with 'if){p, kTg) = H^.c^(p, the aperture-plane wave-front in a 
particular direction in the field and (3i,kTs) the DM-produced correction phase in 

that direction, the mean square residual piston-removed wave-front error over the telescope 
aperture is al = ||efc||^ 2 . 2 (o) ~ where hT is a positive-dehnite weighting matrix that 

removes the piston contribution over the telescope aperture and L 2 is the Euclidean 
norm over the aperture D. 

The integration, reading, processing and correction in the AO loop are done in discrete¬ 
time leading to intrinsic lags which take customarily several milliseconds. We will expand 
the discussion of time delays for synchronous and asynchronous systems in §2.El 

Using the definitions of integrated phase and commands in §2.AI we can now ascertain 
the discrete-time minimum residual phase variance (MV) criterion 

.. M-l 

J(u) = lim — al 

M^oo M ^ 

k=0 
^ M-l 

= \\'^k+iip,(3) ( 6 ) 

k=0 

which neglects any WE dynamics during the integration. The error that is neglected here has 
been dubbed insurmountable error due to the use of averaged variables instead of continuous 
ones and is negligible for the most common systems operating at high frame-rates. 

Provided the phase vector is known the solution to Eq. ([6]) is straightforwardly 
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shown to be a least-squares fit onto the DM influence functions 9|, [15 1 


Ufc(p,/3) = 


(7) 


The latter would have been likewise found from the degenerate control Riccati equation, 
whose discussion we avert here but can be found for instance in QQ. 


We now deal with the general case where V’fc+i is not known in Eq. ([7]) but needs be 
estimated from linearly related noisy measurements. 

The separation principle 0 provides immediately the solution 


Ufc(p,/3) = 

= (N''WN)-Ww5»fc+i|,(p,/3), 


( 8 ) 


where /3) = '0fc+i|5,(c.)(p,/3) = /3)|‘5fc(p, a)) is the conditional 

expectation of '0fc_,_i(/3) in the science directions with respect to the sequence of all 
measurements from (in general non-coincidental) GS directions available up to t = kTg, 
i.e. Sk = jsn, • • • , sM(p, o:) jl8[ |. This conditional expectation is readily available from a 


Kalman filter 


, id, 17|. It relies on the definition of the first-differences state-space model 


in Eq. 


whose construction is shown further below. 


2.C.I. near-Markovian time-progression model of the Ist-order 


The use of auto-regressive models has been given wide attention using the Zernike 


polynomials’ expansion basis set 


4, lid] cis well as with Fourier modes for which a complex 


model encodes perfectly spatial shifts and thus frozen-flow 


20 |. 


We restrict our attention to the near-Markovian first-order time-evolution model 
las delivered best overall performance albeit with increased computational 

, 13 , 


2l| which 


complexity 


22| (spatial dependence omitted) 

V’i:+a(“) = AV>t(Q:) + I'k(a) 


(9) 


with A = r/Tg is a delay in units of sample step Tg, where the transition matrix minimizes 
the quadratic criterion [l|] 


A = argmin{|||V’fc+A(a) - AVfc(a) 11^2(0)) 


( 10 ) 
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( 11 ) 


The solution to Eq. (fTOj) is 

A = i'lpk+A'iPl) (a) 

Note that we do not remove piston in this equation although we could following formulae in 
23|. 

By denoting '4>k{pi kTg) and substituting in Eq. flTT]) we get 

A = ( 12 ) 

= H„ {v(p - ^t)v(pY) {v(p)v(pyy^ H*. 

where 'i/’k+Ai.Py = ^otV^iP ~ ^^kTs) = Hq,'(^(p, kTg) which translates the spatio-angular 
nature of the computation we are accomplishing since a time-translation is factored in as an 
altitude-dependent angular shift using the frozen-flow hypothesis. Figured] shows pictorially 
the statements made. 

Assuming stationarity, the state excitation noise covariance matrix is found for a first 
order time-evolution model from the covariance equality (implicit indices dropped out) 

, (13) 

since (V’fc+iV’I+i) = = (V’V’"'')- The excitation noise covariance matrix is therefore 

= ('0'*/’"'') “ A a”’” (14) 

The model driving noise covariance matrix is a key element of the KF design. 


2.C.2. Measurement model 


The presentation that now takes place restricts itself to the case of the Shack-Hartmann 
(SH) WES which provides the averaged phase gradient over a regular grid of sub-apertures 
in the pupil-plane j24| |. 

We use a geometrical linear model which we assume valid for the regime of operation it 
will be submitted to; the model is thus 


Sfc(a) = GV>fc(Q:)+ 77 ^( 0 :) (15) 

where represents photon an detector shot noise taken to be additive Gaussian noise of 
know covariance and G is a matrix whose entries are computed from the discretization of 








Fig. 1. Spatio-angular computation scheme. Pupil displacements are highly exaggerated for 
viewing. The spatio-angular correlations are computed with the pupil located at p — vr where r 
is time in seconds. 


the wave-front gradients in the GS directions. It uses a 3x3 regular stencil to compute 
the averaged phase gradient with Simpson weights for the fully-illuminated sub-apertures 


Wx = 

a 


1 '^y ~ 


-w^ 


(16) 


-1/4 0 1/4 
- 1/2 0 1/2 
-1/4 0 1/4 

The stencils for partially illuminated sub-apertures at the outer and inner edge of the 
telescope pupil are found by expanding the phase on a set of bilinear splines and computing 
the average gradient over each of the four quadrant patches. The derivation of this result is 
out of scope but can be found at request in [25 1. 

This sampling and respective estimation of wave-fronts at twice the SH-WFS’s Nyquist 
frequency permits the rejection of some aliasing that affects the original measurements as 
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outlined in 


19| . An enhanced model with coloured noise is shown in 261. 


Note that in our implementation we will not make explicit use of the ray tracing H in Eq. 
(IT^ since the layered phase need not be estimated - no multi-conjugation step ensues. Only 
the pupil-plane phase is of interest which opens up to considerable modelling simplifications 
using a spatio-angular frame-work. 

This is a fundamental simplihcation with respect to previous tomographic LQG 
implementations - which in the remainder is called E-LQG to account for the explicit 
estimation step. We avoid estimating the layered wave-front as is customarily done, focusing 
in estimating the integrated pupil-plane wave-fronts in the GS directions along their lines of 
sight. The angular anisoplanatic estimate is made on a subsequent step by invoking linearity, 


stationarity and Gaussianity as in the static case 


Q. 


2.D. Minimal state-space representation: compactness in Open-Loop AO 


The combination of enhanced spatio-angular first-order predictive models with single¬ 
conjugate correction allows for some modelling simplifications that in turn will make for 
smaller problem size and less stringent real-time and off-line computational requirements. 
Under these working assnmptions, the LQG model admits - unlike previons 


implementations jd, l22| - a minimal state representation with phase taken at one single 
time instant. To do so one assnmes the measurements are available at the end of the WES 
integration step, regardless of the actual temporal delay in the system 9j: one affects thus 
the whole delay to the commands instead of the measurements. 

On top of this, we can further use a state vector with the pupil-plane phase integrated 
over directions ck G [1, • • • , N^] at instant k to provide the Spatio-Angular LQG formulation. 

Selecting thus = 'ipkict) 




^ 'ipicxi) ^ 




(17) 


where the state is a concatenation of phase vector in the guide-star directions, one defines 
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the state space terms for Eq. ([9]) 


A 

B 

V 


A 

0 

I 

C 

V 

0 


G 

0 

0 


( 18 ) 


The implementation of the KF involves a real-time state update and prediction equations 
which in the SA case are 


+ "Hoo (sfc(Q:) - 

^k+i\kif^) = ^^k\kioc) 


(19a) 
(19b) 

where 'Hoc is the asymptotic Kalman gain computed from the solution of an estimation 
Riccati equation j9|. The use of the asymptotic value is justified, like elsewhere 271 for one 


seeks long-exposure performance therefore employing the steady-state gain with no loss of 
performance. Whenever implicit we drop the notation (a). 

Under the choice of state made earlier, the L 2 criterion in Eq. (|6]) assumes the equivalent 
form 


J(u) = lim 

M^oo 



1 Xfc 1 

r( 

1 1 

' 1 


A^A -A'^N 
-N'^A N'^N 



( 20 ) 


from which Q = A^A > 0, R = N'^N > 0 and S = -A'^N; this construction ensures 
a stable controller since A is a positive dehnite stability matrix with \eig{A)\ < 1 and 
Q - SRS'^ > 0 and R a full-rank matrix. 


2.E. Time-delayed correction: synchronous and asynchronous commands 

From phase estimates future disturbance can be computed using the predictive 

models above for either the synchronous and asynchronous correction cases. The latter is 
dealt with straightforwardly since now the commands are computed for time step k + A, for 
fractional delays either integer or fractional multiples of the integration time Tg. We further 
assume that due to discrete nature os signals and sequential operations t = Tg + S leading 
to A = r/Tg = 1 -|- S/Tg, i.e. one frame delay plus pure system lag. 

Both cases are shown in Fig. [2] providing the temporal sequence of operations. 
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(k-2)T (k-1)T kT (k+1)T 



S 



Fig. 2. Temporal diagrams. Top: 6 < Tg] Bottom; 6 > Tg. 

The controller is applied in real-time by computing, at iteration k for either case S > Tg 
or 6 < Tg, 

V’fc+i|fc(«) = AV^fc|fc_i(Q:) + AKoo (sk{ct) - (21a) 

Ufc+5/T,(A) = 

= Ts^k+i\k{f^) (21b) 

where estimates the lead wave-front 6 seconds ahead in time to overcome any 
pure lag delay (Sms for the case of Raven) and multiplication by the anisoplanatic 
hlter ^ provides the pupil-plane wave-front 

estimate in the science direction of interest. This two-step scheme is equivalent to a single- 
step prediction on account of the conditional expectation properties, namely 

E{$|^ = s} = E{<h|E{F|R = s}} (22) 

The DM projection is made in parallel for all the science directions by multiplying by 
the generalized inverse of the DM influence matrix in Eq. ([7]). 
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3. Static and predictive Spatio-Angular wave-front estimation 


The state-space framework offers a suitable means to derive static reconstructors against 
which we will compare the LQG formulation. The presentation that follows is a brief review 


of the standard minimisation criteria in AO that can be found elsewhere 


Q, Q, y, 


28|. 


Here we derive it seamlessly from choosing a specihc combination of matrices in the LQG 
controller to set the context for the simulation results that ensue. 


3.A. Static case 

t is a well established result that the static minimum mean square error MMSE solution 

Q 


= (V’pV’I) (G G^ + (W)) So (24) 

can be derived from the LQG solution by taking A = 0 and therefore the conditional 
expectation '4’k\k-i = 0 in Eq. (I2T]) . In such case the solution of the Riccati equation is 
Soo = Sy, and the Kalman gain boils down to 

W'oo = QT (G + {W))~' (26) 

which provides the same resulting estimate phase vector as in Eq. ([23]). 


3.B. Predictive case 


To reduce the lag error, spatio-temporal prediction may be built-in the reconstructor. In 


such case one computes instead 


29| 




(26) 


where ('0fc+A(/3)®I(‘^)) = (V’(P “ vr,/3 )s"''(q:)) is the covariance matrix of pupil-plane 
phase between the science and guide-star directions at two time instants separated by a 
lag r. 
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Equation (126|) can be split into a spatial estimate followed by a temporal prediction step, 
being written in the two-step form 


^fc+A(/3) = AAV’fc(/3) 

(27) 


(28) 


Both this two-step approach or the all-at-once are equivalent as pointed out before due to 
conditional expectation properties. It remains true however many steps are made as long as 
the starting point is a pupil-plane phase estimate. 


4. Computational aspects 


4.A. Off-line computation of Spatio-Angular covariance matrices 


The point-wise covariance matrix computation involves only spatial distances from which 
the spatio-angular covariances matrices are obtained analytically sampling the covariance 
function for von-Karman turbulence 


^ 25/67r8/3 


^1 

5 


Lo ) 


5/6 


K. 


5/6 


27rp 


(29) 


with Lo the outer scale of turbulence, ro Fried’s parameter, T the ’gamma’ function and 
hnally a modihed Bessel function of the third order. Parameter p = \p\ is the distance 
between two phase points in the bi-dimensional plane p = 

For the multiple-layer case, the pupil-plane covariance functions are found from the 
weighted sum over the layers where the weighting is given by the relative strength of layer 
I (also known as the fractional ro). 

The spatio-angular covariance matrices are found by sampling and stacking the analytical 
covariance function in Eq. fl2^ layer-by-layer (since layers are considered independent) 
between points tracing to the WFS sub-apertures’ corners and mid-points whose coordinates 
in any layer axe\{px^^^+6x^hl), {Pyi,n^^ynhi) and the separation vector between sub-aperture 
raster index i on WFS n and sub-aperture j on WFS m at layer I is, 

Pi,j,l ~ iPi,n ~ Pj,m) T ^ 1 ( 6 ^ — Pm)i (30) 

where the hrst term of the right hand side of the equation is the separation between the 
sub-apertures in the pupil and the second term is the global WFS separation vector back- 
projected in the atmosphere. Figure [3] provides a schematic view of that procedure. 
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Fig. 3. Spatio-angular computation scheme. Each column in the covariance matrix is obtained by 
sampling the phase covariance function with distances corresponding to the length of the arrows 
as shown. To comply with the 3x3 stencil in the gradient operator, the corners and the mid-points 
are sampled. 

4.B. Real-time computational burden: scaling laws 

Regarding computational complexity, the SA-LQG formulation offers a reduced number of 
real-time operations compared to the E-LQG formulation. Furthermore it provides less than 
one order of magnitude increase compared to MMSE (static or predictive) reconstructors 
for a ELT-sized case with a multiplexing factor of 20. 

Take Np to be the number of DM or likewise science directions; iVact the number of 
valid actuators per DM; the number of GS directions; ~ 2iVact the number of slope 
measurements per WES; = 4iVact is an under-estimate of the number of phase points 
estimated per reconstructed layer and N\ the total number of layers. Table [U provides the 
approximate number of real-time multiply-and-accumulate (MAG) operations. If we take the 
case of Raven with = 2, Alact = 97, = 3, Ng = 144 we get an increased computational 

complexity of about a factor 15 for the SA-LQG and 30 for the E-LQG controller with 9 
layers. 
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Table 1. Approximate real-time computational complexity. 


Op. / Algo. 

SA-LQG 

E-LQG 

MMSE 


9AsA„ 

9NsNo, + 4A^Ai 


ARoo 

A„A^A,A„ 

AiA^AsA„ 

A/3AactAQ,As 

-^^k\k 

4A^A„ 

4A<^Ai 



NjsNsctN^ 

A^Aact A^ + 4AiA^ 



As an example Fig. 0] shows the scaling complexity for ELT-sized systems (e.g. E- 
ELT MOSAIC [? ], TMT’s IRMOS j^) as a function of the number of science targets for 
^act = 3200 (from a 64x64 DM), Ns = = dA'act,-^! = 10 and Aq, = 6 . For a 

maximum multiplex capacity of 20 science targets both the SA-LQG and the E-LQG are 
only a factor 2 and 3 more demanding than the MMSE, i.e. roughly the same complexity. 
We have only considered the MMSE implemented as a full VMM since there is growing 
evidence that this option maps well into GPUs and should therefore be preferred to iterative 


implementations 


30| . For the time being we exclude from this analysis other considerations 


such as memory bandwidth and pipelinability. 


5. Sample numerical models and simulations 


An end-to-end Monte Carlo wave-optics simulation of the Raven tomograp 


lie system in a 


zonal basis has been set up using the object-oriented OOMAO simulator 


3l| to examine 


the performance of the SA-LQG algorithm in relation to simpler static and predictive 
reconstructors. 


5.A. Model Parameters and Validation 

The £gure-of-merit we use to establish the quality of the wavefront correction is the 
ensquared energy (EE); this is computed as the ratio of the amount of light falling on a 
science camera within a given angle projected on-sky to the total amount of light reaching 
the detector. Strehl ratios are also computed for completeness. 

Table |2] contains the main simulation parameters, selected to reflect the physical 
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Fig. 4. Number of Giga multiply-and-accumulate (GMAC) real-time operations required by the 
LQG when compared to the static/predictive MMSE and the Explicit LQG. 3200 actuators, 6 
WES with 2x3200 slopes per WES, phase estimated over 2x number of actuators per DM, 10 layers 


properties available in the Calibration Unit of our Raven system test-bench 32|]. These 
include a three layer atmosphere located at ground, 5.5 and 11km altitudes, an turbulence 
coherence length vq of 19 cm and an atmospheric outer scale Lq of 40 m. 

Model behavior was verihed by collecting a simulated long exposure science image using 
no AO correction mode. The theoretical estimate of the Full Width at Half Max (FWHM) 
of the long exposure PSF for the baseline atmosphere in J band (1.2 pm) can be computed 
for a hnite outer scale using the expression developed in |33[ for the ratio, Lq/vq > 20. 
With ro = 0.156cm at 0.5pm and Lq = 30m (chosen from median conditions on Manna 
Kea), this computation yields an expected FWHM of 0.3891 arc seconds in J-band. The 
FWHM of a simulated long exposure science image was computed to be 0.3865 arc seconds. 
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confirming good agreement between theory and simulation. Note that this test was carried 
during initial validation of the simulation tools before parameters such as ro were hnalized. 

5.B. Results 

Previous work [l| examined the trade-off between lower WFS frame-rates which reduces 
spatial errors due to noise propagation since measurement’s SNR is higher but increases 
temporal lag error. Correction quality was shown to increase with WFS integration time to 
a point and then begin to decrease as the lag error exceeded the SNR gain. 

Table [3] summarizes the peak performances over the frame-rates for varying GS 
magnitudes. 

These tests show that, as in the modal case, better peak performance can be achieved 
using temporal prediction and this peak occurs at lower frame-rates than the static algorithm 
for each magnitude. Providing reduced noise-propagation the SA-LQG algorithm shows that 
the best performance occurs at approximately the same WFS frame-rate as in the predictive 
case, but that overall performance is signihcantly improved, indicating a reduction in spatial 
error. The data is also plotted in Fig. |5l here it can be seen that by using temporal 
prediction equivalent performance can be achieved with higher magnitude GSs when using 
a lower frame-rate. The LQG can deliver performance equivalent to that of the predictive 
algorithm using GSs which are one magnitude dimmer and equivalent to the static algorithm 
using GSs which are greater than two magnitudes dimmer. 

5.C. SA-LQG V. E-LQG 

In appendix 1^ it is shown that the SA-LQG can be deduced from the E-LQG (under certain 
conditions) in a stationary regime. 

The geometrical nature of the problem, however, is such that we can retrieve the pupil- 
plane wave-front from its layered counterpart through ray-tracing but not the opposite. This 
is patent in the noninvertibility of H which imposes a spatial truncation on the wave-fronts 
over the footprints in the GS direction otherwise represented on larger meta-pupils. 
Simulations have shown that 

(llV’fc(«) - > (||V’fc(«) - (31) 
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Maximum Strehl-ratio for a given magnitude and frame-rate 



Maximum Ensquared-Energy for a given magnitude and frame-rate 



Fig. 5. Static reconstructor vs Stand-alone Prediction and LQG; Strehl ratios (top) and ensquared 
energy (bottom) as a function of NGS magnitude. Left axis shows best achieved performance with 
a static reconstructor compared to the predictive and LQG reconstructors; right axis shows the 
OL-WFS sample rate at which that performance is obtained. 


19 



























i.e., the averaged residual phase variance in the GS directions obtained with the SA-LQG is 
higher than the one obtained with the E-LQG. This difference then translates to the science 
when using the anisoplanatic hlter as follows 

{UtW - ( 32 ) 

The anisoplanatic filter per se adds no further error to the estimation since 

(ll*(/3) - = {\\M0) - > (33) 

A modal decomposition of the errors using Zernike polynomials showed that the low-order 
modes such as tilts, focus and astigmatism are more affected. This is also in agreement with 
the little impact in EE observed which would suffer more from poorer correction of high 
order modes. 

For the Raven case, the differences observed in the science direction are of the percent 

level. 


6. Raven Laboratory Test Results 


Raven’s physical design and system capabilities are summarized in j3|, l3^, l3^; it is equipped 
with a calibration unit (CU) which also serves as a telescope simulator and turbulence 
generator. A broad spectrum source illuminates an array of pinholes, any three of which 
can be picked off to create an NGS asterism, and a series of filters which allow incremental 
changes in the magnitudes of all asterism stars at once. A cross section of measurements 
across NGS WES frame rates were taken for three different hlters to show the improvement 
in performance using prediction and SA-LQG control over the static reconstructor. The 
laboratory results also validate our simulation results which indicate that, using the SA- 
LQG reconstructor, the system can produce performance equivalent to that obtained using 
the static reconstructor when GSs are 2 magnitudes fainter. A plot of the results of these 
measurements is shown in Fig. El The diameter of the asterism used is approximately 2 
arminutes. The CU lamp does not illuminate all pinholes with uniform intensity, therefore 
the data points labeled Mag 13 were generated using GSs with magnitudes [12.6,13.25, 
13.25]; those labeled Mag 15 are actually [14.3, 15, 15], and those labeled Mag 17 are [16.5, 
17.2, 17.2]. 


20 






50 

45 

40 

35 


O) 30 

a> 

c 

lU 

^ 25 

a> 

3 20 

O’ 

v> 

c 

UJ 15 


Ensquared Energy 



- ■© - 
- ^ - 
-«- 
, , o 


LQG Mag 13 

- Prediction Mag 13 

- Static Mag 13 
LQG Mag 15 
Prediction Mag 15 
Static Mag 15 
LQG Mag 17 
Prediction Mag 17 
No correction 


<D 




200 180 160 140 120 100 80 

WFS Camera Frequency [Hz] 


60 


40 


20 


Fig. 6. Measured ensquared energy in a 140 milliarcsecond slit of science images taken while running 
the NGS WFSs at varying frame rates and execnting static, predictive and LQG algorithms at three 
different magnitude settings. The static reconstructor results are not shown for the faintest stars 
as the performance is not better than no AO correction. 

The performance achieved in lab tests is lower than that predicted by simulations in all 
reconstructor cases. However the trend of increasing peak performance with each increase in 
reconstructor complexity is reflected in both the simulation data and the measured data. The 
overall decrease in performance can be attributed to multiple sources; these include imperfect 
calibration, which has a significant impact on OL systems, underestimation of noise sources 
in the simulation compared to reality, effects of the rotated WFSs, DM fitting, OL go-to 
errors, and non-common path aberrations (NCPA) between the OL and science paths, as 
well as between the CL-WFSs and the science camera. A higher amount of total system lag 
than anticipated may also be reducing the performance of all three types of reconstructors. 
Albeit, the main simulation hndings in improved performance could be reproduced with 
the test-bench. We are currently investigating ways of enhancing the calibration and hne 
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knowledge of system parameters in order to optimize overall system performance and gain 
understanding over unmodelled system features. 



arcsec. 


Fig. 7. Comparison of Science camera images for SA-LQG, Predictive and Static tomographic 
reconstruction algorithms. Peak frequency PSF shown for each reconstructor using magnitude 15 
NGSs (see Fig. [6l) . 


The science images in Fig. 0 represent the best performance achieved for each 
reconstruction method for a fixed magnitude; they clearly show that the EE is increased 
and the spot image becomes smaller for both predictive and LQG algorithms over the static 
reconstructor. 


7. Summary 

We have outlined the dynamical Strehl-optimal LQG controller providing minimum residual 
phase variance tailored to MOAO systems. The latter does not require explicit knowledge 
of the tridimensional wave-front prohle allowing for a simplified forward model. Our 
implementation uses anisoplanatic filters to estimate pupil-plane wave-front profiles from 
any other directions based on knowledge of the turbulent prohle. One such approach can 
also be useful for laser tomography AO. 

We have included temporal prediction by assuming frozen-how, in which case time 
evolution can be converted to spatial translations. We have proposed a minimal space-state 
representation with a single temporal instance providing reduced computational complexity. 
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By allowing the controller to run slower while performing temporal prediction, we can 
increase the SNR in the wave-front measurements by collecting more light in each sensor’s 
sub-apertures. We have shown that we could gain two magnitudes for a same performance 
obtained with a MMSE static reconstructor, leading to a sky-coverage improvement of 
roughly a factor 5. 
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Appendix A: Relationships between SA and Explicit LQG 

We set forth to show that the SA-LQG formulation provided in this paper shares important 
features with the more standard “explicit” E-LQG formulation which involves the estimation 
of phase in the layers. 

Let the E-LQG state update equation 

^k+i\k = ^E^k\k-i + (Al) 

where upper/subscript ’E’ indicates matrices for the E-LQG case. 

For the near-Markovian order 1 time-progression model assumed earlier, the driving noise 
statistics are such that — AgSj^Ag. Multiplying out by the ray-tracing operator 

to the left and its transpose to the right equates to — AS.^ A'*' by plausibly taking 

A = HAbH^ (A2) 

and assuming to be computable and furthermore that HHl = I. This assumption 
however may not hold since H is non-squared and its structure gives rise to a rank deficient 
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matrix. We thus expect (and have evidence for) some differences using the SA-LQG with 
respect to the E-LQG formulation. 

We now add the following 

HAe = H(HtAH) 

= AH 

and 

= HAeAI^ 

= ahai^ 

= AAloo 

Multiplying the state update equation on both sides by H to the left equates to 

+ H£^(st - GHvi|t_i) (AS) 

= H(A, - + HOi (A9) 

With the above identities one gets = AAlooG'j/’fcifc-i) yielding 

^fc+i|fc = A(I - AlooG)^fc|fc_i + AlooSfc (AlO) 

which is the SA formulation. 

We are left with the derivation of equivalence of Kalman gains and respective Riccati 
solutions, both computed off-line. The Riccati equation writes 

= AeS^AI + - AeS^H^G^ 

(GSo^G^ + S^) GHS^Ae'^ (All) 

Multiplying by H on the left and H"'' on the right one gets immediately 

5jcxd = ASooA"'' -|- (A12) 

- ASooG'^ (GSooG'^ + S^) GSooA'" (A13) 

where we take Sqo = with the state noise covariance matrix in the SA case obtained 

from the E-LQG using the dehnitions above. In practice this assumption does not hold in 
part due to the finite number of iterations used to compute the asymptotic Riccati estimation 


(A3) 

(A4) 

(AS) 

(A6) 

(A7) 
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error covariance matrix Soo and its E-LQG counterpart, in part since the geometry of the 
problem imposes a spatial truncation on baselines over which covariances are computed 
unfavourable to the SA-LQG. 

The Kalman gain can then be obtained as 

(A14) 

A4oo = So,G'"(GSooG^ + S,)"' (A15) 

which is the straight SA formulation. 

Further to the results shown in Sect. 15.G1 by using the Kalman gain for the SA-LQG 
from Eq. flA14p does not change the performances reported. 






Table 2. Raven Baseline Configuration Parameters. 


Telescope 


D 

8 m 

Atmosphere 


ro 

19 cm 

Lq 

40 m 

zenith angle 

0 deg 

Fractional tq 

[0.596; 0.224; 0.180] 

Altitudes 

[0, 5.5, 11] km 

wind speeds 

[5.68; 6; 17] m/s 

wind direction 

[90; 180;180] deg 

Wavefront Sensor 


Nngs 

3 

Order 

10x10 

Spix 

0.4 arcsec 

RON 

0.2 e" 

NGS radii 

45 arcsec 

^pix 

12 

f sample 

30-200 Hz 

^WFS 

0.64 /xm 

Centroiding algorithm 

thresholded Centre-of-Gravity 

DM 


Ndm 

2 

Order 

11x11 

stroke 

infinite 

influence 

cubic 

AO loop 


pure delay 

Tag = 3ms 

Evaluation Wavelength 


^evl 

1.65 fim 
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Table 3. Raven End-to-End simulation results. The optimal performance ( % ensquared energy) for 
each GS magnitude is shown for the zonal static SA, and compared to stand-alone SA Prediction 
and the full SA LQG algorithm. 


GS mags 

static SA 

EE lag Strehl lag 

SA Prediction 

EE lag Strehl lag 

EE 

3A LQG 

lag Strehl lag 

13.5 

47.99 

6 

28.55 

13 

49.88 

23 

31.70 

15 

52.27 

20 

33.54 

15 

14 

47.48 

10 

28.05 

14 

50.01 

25 

32.32 

22 

52.66 

22 

34.60 

22 

14.5 

46.79 

12 

27.98 

15 

49.65 

24 

31.72 

24 

52.48 

24 

34.02 

23 

15 

45.73 

13 

25.64 

15 

49.07 

25 

30.66 

25 

52.06 

25 

33.12 

25 

15.5 

44.36 

15 

23.41 

18 

48.41 

26 

29.11 

25 

51.40 

25 

31.53 

25 

16 

42.57 

16 

21.52 

22 

47.31 

30 

26.50 

25 

50.27 

26 

29.25 

27 

16.5 

40.43 

20 

18.90 

25 

46.73 

41 

24.33 

40 

48.98 

35 

26.72 

35 

17 

38.28 

25 

15.57 

25 

45.20 

46 

23.29 

47 

47.48 

45 

24.95 

45 


30 








